Theory of Radio Frequency Spectroscopy of Polarized Fermi Gases 
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We present two exact results for singular features in the radio frequency intensity I(ui) for ultracold 
Fermi gases. First, in the absence of final state interactions, I(u)) has a universal high frequency 
tail Clj- 3/2 for all many-body states, where C is Tan's contact. Second, in a normal Fermi liquid 
at T = 0, I(u>) has a jump discontinuity of Z/(l — m/m*), where Z is the quasiparticle weight and 
m* /m the mass renormalization. We then describe various approximations for 7(w) in polarized 
normal gases. We show why an approximation that is exact in the nj. = limit, fails qualitatively 
for ni > 0: there is no universal tail and sum rules are violated. The simple ladder approximation 
is qualitatively correct for very small nj,, but not quantitatively. 



There has been intense experimental activity on char- 
acterizing, various states of matter in ultracold atomic 
gases [l], 0] ■ This will become ever more important with 
the possibility of new and exotic states being realized in 
these systems. An important tool in these studies is radio 
frequency (RF) spectroscopy where an RF pulse is used 
to transfer atoms from one hyperfine level to another. 
The RF signal [H, BJ has turned out to be much harder to 
interpret than initially thought because of complications 
of strong final state interactions and the inhomogeneity of 
trapped gases. Recently it has become possible to elimi- 
nate these problems by choice of suitable hyperfine levels 
(in 6 Li) and by tomographic techniques that focus on spe- 
cific regions of the gas. The most detailed experimental 
results are available for polarized Fermi gases 0, . 

Motivated by these experiments, we first describe two 
exact results for the singular features in RF spectra of 
a two-component Fermi gas with arbitrary interactions. 
Our results complement the exact results on sum rules 
@j Q ■ We work in the limit where final state interaction 
effects are negligible, so that we can focus on the non- 
trivial effects of interactions in the many-body state. 

(1) The RF spectrum I(uj) has a universal Cw -3 / 2 tail 
at high frequencies, where C is Tan's contact coefficient 
0, which is independent of spin. This form is valid for 
all phases of Fermi gases: super fluid E l 3 , hig hly unbal- 
anced normal Fermi liquid [y, Iff! fTlTTl3 . If3l | or even a 
balanced Galitskii Fermi Liquidfl^. 

(2) In any normal Fermi liquid state, the RF spectrum 
I(lu) at T = has a jump discontinuity. Its location de- 
pends on the chemical potential /i and its magnitude is 
determined by the combination of Fermi liquid parame- 
ters Z/(l— m/m*), where Z is the quasiparticle weight 
and m* the effective mass. 

These exact results are important not only in inter- 
preting experiments, but also in understanding various 
approximation schemes 0, [TH [l6[ , which are necessar- 
ily required to calculate the RF lineshape I(u>) for a 
strongly interacting gas. In the second part of our pa- 
per we critically analyze diagrammatic approximations 
for the highly unbalanced normal Fermi liquid. 

(3) We show that a simple self-consistent approximation, 
motivated by the fact that it is essentially exact in the 



nj_ = limit [l3T |. has serious qualitative problems for 
non-zero n^: The minority spins do not exhibit the uni- 
versal tail leading to sum rule violations and majority 
spins are completely unaffected. 

(4) A simple ladder approximation, on the other hand, 
correctly exhibits all of the qualitative features expected 
on general grounds for rij > 0, however there are quan- 
titative inaccuracies and the approximation breaks down 
for ni > 0.05. 

Formalism: Consider a Fermi gas with three hyper- 
fine states which we label as |, j and "e" (for excited 
or empty). The number density in levels a =T,| is n a 
with corresponding (non-interacting) Fermi energies tFa- 
The | and J. fcrmions interact with an s-wave scattering 
length ay, i = a. The e-level is located at energy AE a 
from the bottom of the a bands, and is empty (n e = 0). 
We assume that fermions in the e state do not interact 
with those in a levels: a e (T = 0. If such interactions 
are strong, the simple results obtained below are consid- 
erably modified by vertex corrections (l6| . One is then 
dealing with the complications of the probe in addition 
to the many-body system of interest. 

When final state interactions arc negligible linear re- 
sponse theory leads to the simple result: 

I a (ui) = y^ a (k, e k - jjL a -cj)n F (e k - \i a - ui) (I) 



where the RF shift ui = u)rf — AE a , ljrf is the RF 
frequency, e k = k 2 /2m is the bare dispersion, \i a the 
chemical potential, and h = 1. rii?(e) is the Fermi func- 
tion and the single particle spectral function A a (\s.,ui) = 
— ImG CT (k, uj + i0 + )/7r includes all many-body rcnormal- 
izations due to interactions between j and J, fcrmions. 

Sum Rules and Large-w behavior: The exact sum 
rules @, 111 for the zcroth {I = 0) and first (1 = 1) mo- 
ments of the RF intensity J diuui I a (oj), are valid for all 
values of a and a e ^ a . It might seem that the first moment 
sum rule (clock shift), which diverges as a ej(T — > 0, can 
be of no use when final state interactions are negligible. 
However, we find that this divergence is actually related 
to a universal high frequency tail in I a (uj). 

We rewrite (I) as -^(w) = 
J2 k J dnA a (k,Q)n F (n)S(n - e k + n„ + w). This 
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immediately leads to zeroth moment sum rule |17j | 
/ duil a (uj) = N a , using / d£lA a (k, fi)np(fi) = n a (k) 
and n CT (k) = N a is the number of a fcrmions. 

We next analyze l a (io — ► oo). The delta-function 5(0— 
£k + Mff + w ) then contributes in one of two ways: cither 
(a) is large negative with eu small, or (b) small but 
£k large. In case (a), however, the spectral function A a 
vanishes for small k and fl — > — oo. Thus only case (b) 
survives and we find I a (uj — > oo) ~ ^ k 7i CT (k)(5(ek — £*>)■ 
Using Tan's result [§] n^k) ~ C/k 4 for k ^> kp we thus 
find that 



(2) 



where C is the contact. We emphasize that the form 
of this result is independent of the phase (normal or su- 
pcrfluid) of the Fermi gas, though the value of C does 
depend on the phase. (This tail is absent only for the 
noninteracting gas for which C = 0.) Note that this high 
frequency tail arises from short-distance physics in any 
Fermi gas, and is crucial for enforcing the divergent clock 
shift for a e!(T = 0. 

Fermi liquid singularity: In the study of many- 
body systems, various phases are often directly identified 
by characteristic low-energy singularities in measurable 
quantities, such as the the discontinuity at kp in the mo- 
mentum distribution of a Fermi liquid, or the square root 
singularity in the density of states of a s-wave supercon- 
ductor at T = 0. Here we ask if any such singularity 
exists in the RF signal. Given the k-sum and the kine- 
matics in eq. (1), we see that there is no characteristic 
singularity in the paired superfluid state. However, as we 
show next, there is a singular signature for normal Fermi 
liquids at T = 0. 

In the remainder of this paper we focus on the normal 
(i.e., non-superfluid) ground state of the highly polarized 
Fermi gas. Thus our results are relevant, e.g., to the 
unitary gas which has been predicted to be a normal 
Fermi liquid for x = n\ /rtt < 0.4, based on quantum 
Monte Carlo simulations [l^]. (Our general results a pply 
equally well to the dilute repulsive gas of Galitskii [14| . 
which is yet to be realized in the laboratory.) 

For a Landau Fermi liquid the spectral function is 

A(k,v) ~ ZS(w- k F (k~ k F )/m*) + A inc (k,u) (3) 

close to the Fermi surface (k ~ kp,u> ~ 0). The sub- 
script a is dropped for simplicity. The first "coherent" 
term gives the quasiparticle pole in the Green's function 
with quasiparticle weight Z and effective mass m* [is| . 
kp is unshifted from its bare value as required by Lut- 
tinger's theorem [Til ]. The second non-singular term is 
the "incoherent" part of the spectral function. 

The singular contribution to I(u>) is obtained by sub- 
stituting the coherent term in ([3]) into |T]) and using 
np(e) = Q(— e) at T = 0. We convert the k-sum to 
an integral over ek and write the quasiparticle dispersion 
as k F (k — k F )/m* ~ (fc 2 - kp)2m* = (e k - e F )m/m*. 



For m* > m we find a peak which grows like a square 
root in u> and then has discontinuous drop, all of which 
rides on top of top of the smooth contribution from the 
incoherent piece. The location of the discontinuity u>* 
and the size of the jump AI are thus given by fl9j 



Ah = 



ZgN{e Fa ) 
(1 — m/m* ) 



(4) 



where N{eF a ) is the density of states at the Fermi energy. 
Diagrammatic Lineshape Calculations: The 

form of the Fermi surface singularity and the high 
energy tail in the RF intensity have been elucidated 
above on general grounds. Calculating the detailed 
lineshape necessarily requires approximations 

to be made for a strongly interacting Fermi sys- 
tem. Here we describe diagrammatic calculations 
for the highly imbalanccd normal gas, highlighting 
the successes and limitations of two approximation 
schemes. All such calculations sum particle-particle 
(p-p) channel ladder diagrams: r _1 (q, iqi) = m/Ana — 
E k [l/2e k - /T 1 £„ Gj (k + q, ik n + iqe)G L (-k, -ik n )] . 
The rationale for focusing on p-p ladders can be given 
in many different ways. These are the leading diagrams 
for Fermi systems with short-range interactions in the 
low density repulsive Fermi liquid [l4| . in the Nozieres- 
Schmitt-Rink analysis [2(| of the normal state of the 
BCS-BEC crossover, and in the 1/N expansion for the 
attractive Fermi gas [l5| . 

One can analytically obtain closed-form expressions for 
the real and imaginary parts of the retarded r _1 (q, to + 
i0 + ) when G a are the bare Green's functions; details 
are omitted for simplicity [2l[. Next, the self-energies 
E <y (k,ifc n ) = /3 _1 J2 u jT(q,iqi)G^ a (-k + q, -ik n + 
iqe) are calculated [22| using the spectral represen- 
tation for the p-p vertex in terms of ImT(q, to + 
i0+) to obtain E^k,^ + i0+) = E' CT + This 
in turn leads to the spectral functions A a (k, uS) = 
— Im [lu — ek + — T/ a — /it, which form the ba- 

sis for our calculation of n(k) and of the RF spectrum 
using (1). 

(I) Let us first discuss a simple self-consistent approxi- 
mation Q , motivated by an analysis that reproduces the 
essentially exact result [13| of single J, spin (n\ = limit) 
interacting with a Fermi sea of f fermions [12J ■ We will 
show that this scheme has serious qualitative problems 
for > and analyze why this is the case. In this ap- 
proximation, the Green's functions used to calculate T 
and E are the bare G's but with a renormalized chemical 
potential. A self-consistency condition is then imposed 
so that /ij = £fi + E^(&f|, 0; /ijj, where Ej itself de- 
pends on fj,i [l3j . For the single minority spin limit, this 
reproduces the result £tx( n l = 0) = ~ — — 0.6ei?f. 

This approximation for rij > implies the use of a 
negative fx^ in the bare G± used in T and E. As a re- 
sult, one misses all effects of finite uj. occupancy "in- 
side" the calculation. We can then analytically see that 
Imr(q, lu < 0) = which impacts the results as follows. 
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FIG. 1: (color online) RF spectrum of a unitary Fermi gas 
with ni/n-[ = 0.25 calculated within the self consistent ap- 
proximation (I) (see text). The majority (f, red) is a delta 
function with weight nj. The minority (J., blue) spectrum has 
a discontinuity and a shift due to interactions, but no high 
frequency tail. 



FIG. 2: (color online) RF spectra of a unitary Fermi gas with 
n l/ n T = 0.05 within a simple ladder approximation (II) (see 
text). Both minority (J., blue) and majority (|, red) spectra 
exhibit a discontinuity (dashed) and a large-u tail. The blue 
curve for minority spins is 15 x I\.(u)). 



For the minority fermions £"(k, u < 0) = 0, which im- 
plies Af c (k,oj < 0) = and thus uj.(k) = Z l Q(k Fl - k). 
This means that X^k^lW = ZNi < N± and the zeroth 
moment sum rule for Ii(u>) is violated. In addition, in 
the absence of any incoherent spectral weight for u> < 0, 
one also misses both the universal k~ A tail in nj.(k) and 
the w -3 / 2 tail in the RF spectrum (see Fig. [T]). The first 
moment of J'j.(w) is then finite, instead of diverging as it 
should. Further, the majority spins are completely unaf- 
fected by interactions in this approximation since one can 
see analytically that Ej-(k, w) = 0, which is clearly un- 
physical for non-zero ny. The majority (f) RF spectrum 
is thus a delta function, a result that is at odds with all 
available experiments. Clearly this approximation fails to 
provide a reasonable description of RF spectra of highly 
imbalanced gases, despite its success in obtaining reason- 
able numerical estimates for hi . All of the problems here 
arise from the fact that propagators with renormalized 
Hi < are used without taking into the shifts in the J. 
particle dispersion. 

(II) This suggests that it may be physically more sen- 
sible to do the simplest calculation without any attempts 
at partial self-consistency, i.e., evaluate all diagrams with 
bare propagators and bare (.Fa- This leads to equations 
which are identical with the 1 /N approximation [15J with 
N = 1 at the end. Now, in contrast to the previous ap- 
proximation, the p-p vertex ImT has structure even for 
lj < 0, and this leads to n(k) and I(oj) with universal 
tails for both spins. 

Our numerical results for the RF spectra of a highly 
polarized unitary Fermi gas with x = nj,/ n t = 0.05 are 
shown in Fig. O Both majority and minority spectra 
show jump discontinuities and high frequency tails. The 
spin-independent oj~ 3 / 2 behavior tails are observed in 
Fig. [3J Comparing the results of (I) the self-consistent 
Ha approximation in Fig. [T] and (II) the simplest ladder 
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FIG. 3: (color online) High frequency tails of the RF spectra 
of the unitary Fermi gas with ni/n-f — 0.05 shown in Fig. [2] 
Both majority (red, open circles) and minority (blue, filled 
circles) spectra exhibit a lj~ 3 ^ 2 tail. 



approximation in Figs. [2] and [3l there is no doubt that 
the latter provides a far better qualitative description of 
the RF spectrum. 

Despite these qualitative successes, it must be empha- 
sized that the simple ladder approximation (II) is not 
quantitatively accurate insofar as the calculated chemi- 
cal potentials, e.g., Hi = e Fi + ^i(^Ff) 0j e Ff)- I n par- 
ticular we find that in the single spin limit Hl( n l = 0) — 
—0.96fi, as compared with the exactly result of — 0.6eF|- 
Moreover, we have found that the simple ladder approx- 
imation leads to a negative compressibility for x > 0.05 
clearly signaling the limitations of the approximation. 

The prospects for a better diagrammatic approxima- 
tion are unclear, since fully self-consistent calculations 
do not necessarily lead to better answers in strongly in- 
teracting systems [23|. We also note that, while Quan- 
tum Monte Carlo (QMC) calculations have often pro- 
vided valuable quantitative information [l^] for energet- 
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ics, extracting frequency-dependent correlation functions 
from QMC is very difficult in view of two serious issues: 
the fermion sign problem in polarized systems and the 
problem of analytic continuation. 

Comparison with Experiments: Let us begin with 
the universal high frequency tail. A long tail is visible 
in all of the published spectra (Fig. 1 of ref. [f|; Fig. 2 
of ref. Q). It is present for superfluid as well as normal 
state spectra and the same for both spin species, as we 
predict. It would be interesting to know if the signal 
to noise ratio in experiments is sufficient to test the 3/2 
power law [24[ and determine the coefficient C fl7l |. 

The jump discontinuity in the T = RF signal for a 
normal Fermi liquid will be broadened by finite temper- 
ature and by experimental resolution. The best we can 
expect then is to see a peak at, or very close to, the loca- 
tion of the discontinuity. For the minority spins eq. (|4|) 
predicts this to be wT/eFt = 3A /5 — (1 — m/m^x 2 / 3 — 
6Fx/5, using the best QMC result [25[ for This is 
exactly the expression used for the peak position by Schi- 
rotzek et al. [6J|. For the majority spins, with m* ~ m, 
the peak will be at ui^/eF] = 2Aox/5 + Fx 2 /5, which is 
slightly shifted from zero. 

Conclusions: We have derived two exact results for 



singular features in the RF spectra of Fermi gases. The 
high frequency w~ 3 / 2 is a universal feature, independent 
of the nature of the many-body state, when final state in- 
teractions are negligible, and provides an opportunity for 
measuring Tan's contact C . Such a study combined with 
other experimental probes such as photoassociation (2(| 
could provide deeper understanding of how short range 
physics controls the universal properties of strongly inter- 
acting cold gases. Our second exact result on the jump 
discontinuity in the spectrum of a Fermi liquid at T = 
provides a distinguishing feature between a normal and 
superfluid ground state. In the second part of our paper 
we show that, in the absence of a small parameter, it 
is very difficult to obtain reliable results for the detailed 
frequency dependence of the RF spectra - which capture 
both general qualitative features and are quantitatively 
accurate - in strongly interacting quantum gases. Indeed, 
that makes exact results such as sum rules and the singu- 
lar features derived in this paper all the more important 
in interpreting experiments. 

Acknowledgments We acknowledge discussions with 
R. Dicner, W. Ketterle and M. Zwierlein, and support 
from NSF and ARO. VBS thanks DST, India for support 
through a Ramanujan grant. 



[1] 

[2] 

[3] 
[4] 
[5] 
[6] 
[7] 
[8] 

[9] 



[10] 

[11] 
[12] 



[13] 
[14] 



[15] 
[16] 

[17] 



Zwerger, Rev. Mod. Phys. 
Stringari, Rev. Mod. 



I. Bloch, J. Dalibard and W. 
80, 885 (2008). 

S. Giorgini, L. P. Pitaevskii and S. 
Phys. 80, 1215 (2008). 
C. Chin et al, Science 305, 1128 (2004). 
H. Schunk et al, Science 316, 867 (2007). 
Schirotzek et al, Phys. Rev. Lett. 101, 140403 (2008). 
Schirotzek et al, arXiv:cond-mat/0902.3201vl. 
Baym et al, Phys. Rev. Lett. 99, 190407 (2007). 
Punk and W. Zwerger, Phys. Rev. Lett. 99, 170404 
(2007). 
S. Tan, 



C. 
A. 
A. 
G. 
M. 



arXiv:cond-mat /0505200 



Braaten and L. Platter, Phys. Rev 
(2008). 

M.W. Zwierlein 
G. B. Partridge 
S. Pilati and S. 



and 
Lett. 



0508320; E. 
100, 205301 



et al, Science 311, 492 (2006). 
et al, Science 311, 503 (2006). 
Giorgini, Phys. Rev. Lett. 100, 



C. Lobo et al, Phys. Rev. Lett. 97, 



030401 
200403 



Lett. 98, 180402 (2007); 
Phys. Rev. Lett. 101, 



(2008) 
(2006) 

R. Combescot et al, Phys. Rev. 
R. Combescot and S. Giraud, 
050404 (2008). 

A. A. Abrikosov, L. P. Gorkov and I. Dzyaloshinski, 
Methods of Quantum Field Theory in Statistical Physics, 
(Dover, NY, 1963). 

M. Villette et al, Phys. Rev. A 78, 033614 (2008). 

A. Perali, P. Pieri and G.C. Strinati, Phys. Rev. Lett. 

100, 010402 (2008). 

We will use the zeroth moment sum rule N a to fix our 
normalization of I a {uj), since we have set the RF coupling 
constant in (1) to unity. 



[18] The quasiparticle weight and effective mass are related 
to the real self-energy E' via the standard results: Z — 
1/ [1 - dE'/dcj] and m/m* = Z [1 + <9E'/<9e fe ] evaluated 
at the Fermi surface: (k — kF,u = 0). 

[19] If m* ~ m, but Z < 1, the singular contribution to I(u>) 
is a delta-function at u*, and together with a incoherent 
piece which has the universal high frequency tail. 

[20] P. Nozieres and S. Schmitt-Rink, J. Low Temp. Phys. 59, 
195 (1985). 

[21] P. Nikolic and S. Sachdev, Phys. Rev. A 75, 033608 
(2007). 

[22] The iqi Matsubara sum can be done analytically using 
the spectral representation. The q-integrals in E and k- 
integrals in (1) have to be done numerically. 

[23] See, e.g., R. B. Diener, R. Sensarma, and M. Randeria, 
Phys. Rev. A77, 023626 (2008). 

[24] We expect that small a CTi( . = a' will lead to crossover 
from uj- 3/2 to an lj~ 5/2 tail for u > h 2 /ma' 2 . In ref. 
l/(fcp^a') 2 ~ 20. Such a crossover is consistent with the 
exact solution in the Bose limit by C. Chin and P. S. 
Julienne, Phys. Rev. A 71, 012713 (2005). 

[25] We use E= § A T e FT (1 - A x + (m/m* )x 5/s + Fx 2 ), the 
QMC parametrization [ll| for the ground state energy. 
Here x = nj/rif, E^ — 3Ao£ft/5 is the binding energy 
and mj the effective mass of a single J. spin, and F a 
Fermi-liquid correction. We use mj here since m* is in 
general a-dependent. 

[26] G. B. Partridge, et al, Phys. Rev. Lett. 95, 020404 
(2005). 



